Controllable diffusion of cold atoms in a 
harmonically driven and tilted optical lattice: 
Decoherence by spontaneous emission 

Navinder Singh 

Institute of Physics, Bhubaneswar-751005, India 
E-mail: navinderOiopb . res . in 

Abstract. We have studied some transport properties of cold atoms in an accelerated 
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1. Introduction 

The seminal review paper |JJ of Stig Stenholm on the theory of laser cooling gave a 
boost to cold atom physics and optical laser technology [21 [3]. Cold atoms on optical 
lattices constitute clean quantum systems as compared to solid state systems due to very 
less scattering and decohering effects, and provide a paradigmatic model to test some 
novel quantum phenomena predicted half century back for metallic systems (H [5]. The 
experimental advancements has further provided much more deeper understanding in 
fundamental quantum statistical laws, new insights in quantum computing technologies, 
quantum phase transitions (superfluid to Mott-insulating phase) in the regime of strong 
correlations [6], and Bose-Einstein condensation on optical lattices [H [8j [9j [10l ITT] . 

In the present contribution, we consider single atom transport (in one dimension) 
on a tilted optical lattice in the presence of an external AC field. To fix ideas, consider 
the physical situation of an atom interacting with laser field. We have a single atom 
A, with an excited state |e) and a ground state \g) separated by an energy interval 
E e — E g = Tiwa- The atom Hamiltonian is Ha = p 2 /2m + hu^le) (e\ (with E g = 0). 
The atom is subjected to a classical laser field with electric field E(x, t) = eE(x)e~' luJLt , 
where ujl is the laser frequency and e is the polarization vector of the laser. If the 
amplitude of the electric field E(x) is varying slowly in space x compared to the size of 
the atom, the atom- field interaction can be described in dipole approximation, i.e., by 
the coupling = — /iE(x, t), where \i is the atomic dipole moment. We assume that laser 
is far detuned from any optical transition in the atom. In this simplified picture one uses 
perturbation theory and eliminate internal atomic states from the dynamics to obtain 
an effective potential V(x). The atom Hamiltonian thus becomes Ha = P 2 /2m + V(x), 
with V(x) = \Q(x)\ 2 1 'A(ld x — uja)- The term Q(x) = —2E(x){e\[ie\g) is called the Rabi 
frequency which drive the transitions between the two atomic levels. Now consider 
that the atom is in a non-resonant standing light (E(x,t) = 2eE$ cos k^x cosc^t = 
eE [cos(uit — k L x) + cos(uj L t+k L x)], two counter propagating waves) . This constitutes 
'The Basic' optical periodic potential U(x) = Uq cos 2 (/clx), which can be accelerated by 
frequency chirping. Lattice Potential is tilted in the reference frame of the atom (co- 
moving frame), and the atomic motion is governed by quantum mechanics of a particle on 
a periodic lattice (Figure 1). When there is no acceleration (no tilt), an initially localized 
wavepacket will spread through resonant Bloch tunneling and become delocalized. But 
in a tilted lattice (optical lattice being accelerated analogously electron and crystalline 
lattice is in an external constant electric field) atoms can remain localized due to 
suppression of Bloch tunneling. They exhibit novel quantum phenomenon of Bloch 
oscillations due to the repeated Bragg scattering[12j. In the case of electrons in metals 
this corresponds to an induced AC current with applied DC voltage across the sample. 
But in usual practice this purely quantum effect is overshadowed by scattering processes 
and we obtain ohomic DC current. The second interesting effect of potential tilt is that 
the Bloch bands are broken up into Wannier-Stark (WS) Ladders of states [T3l IT4] . The 
level spacing between two nearby levels in the ladder is given by eEd (e is the electron 
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charge, E is the applied constant electric field and d is the lattice spacing, for the 
case of optical lattice, level spacing = Fd = mad, a is the imparted acceleration, and 
d = n/ki is the period of optical potential). As the tilt of potential per lattice spacing 
becomes comparable to well depth, a new interband tunneling process called Landau- 
Zener tunneling becomes important which is a natural extension of the stark effect for 
a single atom. 

Now the above stated phenomena are purely quantum in nature. But the relaxation 
processes are natural. No system is an ideally decoupled system for all space and 
time scales. In case of cold atoms on optical lattices the main relaxation process 
is the spontaneous emission of photons by excited atoms. Since photons has finite 
momentum, it's generation gives a recoil kick to the atom (mechanical effects of light). 
These relaxation processes decoher pure quantum effects. The decay of Bloch oscillations 
and the diffusive spreading of the atoms is shown to cause by spontaneous emission|15]. 
Interestingly, the presence of relaxation processes are actually important for practical 
purposes, like ohomic current across metals(Joule heating effects) and atomic current 
across an optical lattice etc. [16]. 

In the present contribution we consider the effect of external AC drive on the 
quantum transport of cold atoms on tilted optical lattice, and the effect of spontaneous 
emission that cause decoherence, which is essential for diffusive motion. We have 
obtained a tunable diffusion co-efficient and it's nonlinear enhancement with increasing 
drive amplitude and we also report a novel maximum diffusion condition. The analytical 
results obtained by us correctly specialize to the exact results known in the limit of 
zero drive and zero-bias. Also, in addition to cold atoms, the results obtained are 
applicable to experimentally realizable super lattice hetero-structures that support the 
Stark- Wannier (SW) ladder states in the presence of a strong longitudinal electric-field 
bias[H [131 [TU, [21]. As is well known, a strong field normal to the superlattice planes 
can break up the extended Bloch-like band continuum into energetically well resolved 
states localized in the potential wells. The stronger the biasing field the more localized 
are the SW state [331 EE]. 

This paper is organized as follows. In section 2 we introduce the model Hamiltonian 
and the simplified master equation for the present system. In section 2.1 we consider 
the simple case of no acceleration and no external AC drive, and in section 2.2 tilted 
lattice in the presence of an external AC drive. We end with brief discussion of results. 

2. Model Hamiltonian and Simplified Master equation 

We have three interacting physical systems (1) the atom (2) the laser field, and (3) 
the vacuum. Up to this point we have not considered atom-vacuum field coupling, 
which is responsible for spontaneous emission of photons by the excited atom. The 
random recoil kicks (mechanical effects of spontaneously emitted photon momentum) 
cause decoherence of atomic motion, which is the main object under study. Spontaneous 
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BROKEN BLOCH BANDS ( WANNER - STARK LADDERS OF STATES) 




Figure 1. An atom in an accelerated optical periodic potential (in co-moving frame 
of reference it is a tilted or washboard potential) 

emission is characterized by the natural width of the excited state or the radiative 
life-time of the state (=1/7). If we are interested in very short interaction times 
(t <C I/7), we can neglect spontaneous emission, and the evolution of atom-laser 
system is described by Schrodinger equation. But for long interaction times (t ^> I/7), 
due to the presence of several spontaneous emissions, one cannot use Schrodinger 
equation. In this situation one uses system-environment approach i.e., the reduced 
atomic evolution (traced over infinitely many vacuum field (environment) degrees of 
freedom) is then given by a master equation. In the present situation of far detuning 
(5 = — ^0 ^ ^(atomic Rabi frequency) ,the time evolution of the reduced density 
matrix p of the atomic motion along the x direction is given by a simplified master 
equation[15] (detailed theory is given inpH [T8| [19], |2j [20j ) : 

^ = -l[H,p]-l J dup{u)[L{L u p + P L[L U - 2L u pLl] , (1) 

Where 7 = 7o^r(7o is the inverse radiative life time of the excited state), f2 is the 
atomic Rabi frequency, and 5 is static detuning. The first term on the RHS gives the 
unitary evolution, while the second term gives the non-unitary (incoherent) evolution 
causing the initially pure density matrix (p = p 2 at t = 0) to become mixed (p 7^ p 2 for 
t > 0). Here p(u) is the angle distribution of the spontaneously emitted photons, and for 
linearly polarized light, p(u) ~ 1/2. The operator L u is the projection of photon recoil 
operator along the atomic direction of motion. The photon recoil operator represents 
the coupling of the internal atomic dynamics (decay processes) to the external atomic 
motion(center of mass motion). In dipole and rotating- wave approximation, 

L u = cos(k L x)e iukLX \u\ < 1 (2) 

The physical effect of the photon recoil operator is to randomly change the atomic 
quasimomentum. Initial delta function distribution of atomic quasimomentum will be 
smeared over the entire Brillouin zone and cause the decay of Bloch oscillations as 
observed in[15j. Our approach is based on one band tight binding hamiltonians similar 
to|15](new feature added is the AC drive, section 2.2). In matrix element notation, our 
system Hamiltonian and Recoil operator is 

H mn = — — [S m , n +i + ^m,n-i] + Fdm5 mn (3) 
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Figure 2. No acceleration case, all states having same energy. 



L mn = (-l) m e™5 mn (4) 

Here, V(> 0) is the transfer matrix element between nearby states, F is the inertial 
force acting on the atom, d = rc/k L is the period of optical potential, and m is the 
localized Wannier function \m) associated with the Zth well of the optical periodic 
potential. In the following sections we will proceed step by step, starting with un- 
accelerated potential. 



2.1. No acceleration and no AC Drive 



We begin by considering first the simplest case of quantum motion of an atom moving 
on an optical lattice (Figure 2) under a tight-binding one-band Hamiltonian 

/J° = -(\//2)^(|/)(/+l| + |/ + l)(/|) (5) 
i 

where (-V) is the nearest-neighbor transfer matrix element, and the sum is over the 
N sites with N taken to be infinite. The effect of spontaneous emission, namely the 
incoherence, will be introduced through recoil operator defined in equation 4. The 
reduced density matrix p for the particle then obeys the evolution master equation 
(equation 1). In terms of matrix elements, we have 

dpmn _ _iV_, 
dt ~ 2h 
With the initial condition, 



[Pm,n+1 + Pm,n-1 ~ Pm+l,n ~ Pm-l,n] ~ 1 Pmn]}- ~ 0~n 



Pmn(t — 0) — SmoSnQ. 

In Fourier space, with (5 — — jr. 
d 



^ 1 Pm,nt 



-imk\ ink2 



m,n 



with 



We get 



- 7 Pm,ne- tmkl e mk > + 7 £ 5 m , n e-^e mk \ 

m,n m,n 

p(k l7 k 2 ,t) = J2Pm,ne- imkl e ink2 ,5 m , n = ^f e l{m - n)q dq. 



d_ 

dt 



p(k l7 k 2 ,t) = [i(3(cosk 2 - cosfci) - ^f}p(k l7 k 2 ,t) 



i r ~ 

+ 2n J P ( y k ^- ( li k i- < li t ) d< l- 



(6) 
(7) 



(8) 



(9) 



(10) 
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Defining center of mass and relative wave- vectors as p = (k\ + k 2 )/2, u — k\ — k 2 and 
writing p(k±, k 2 , t) = p(p, u, t) we have 

d_ 
dt' 

We define further reduced density matrix by 

1 



7 r 

-p{p,u,t) = [2iPsmpsm(u/2)-~f]p(p,u,t) + — J p(p- q,u,t)dq. 



(11) 



X{u,t) 



2ir 



p(p-q,u,t)dq 



2vr 



p(q,u,t)dq. 



(12) 



As the dimensions of (3 and 7 are (time) 1 , we define the dimensionless parameters as 
t — t/3 and r = ^, and defining 

4>(p, u) = 2isinpsin(-u/2) — T. (13) 

So with this, the evolution equation becomes 



d 

— p(p, u, t) = (f)(p, u)p(p, u, t) + Fx(u, t). 



(14) 



The detailed calculation of mean-squared displacement from the above equation is given 
in Appendix A. We finally obtain 

(AT)) = -^ + fT + ± I e^. (15) 



r 2 



Recalling, r — t(3 and r = 7//?, where (3 = V/h, 7 = 7o§r, with this 

"o 

(^)) = -^ + ^t+^ , 

which reduces to the classical case in the large time (t I/7) limit as (x 2 (t)) ~ ((3 2 / , ~f)t 
giving the diffusion coefficient 



(16) 



D 



(3^_ _V*_ l co L - co 

27 2 l0 h 2 [ n 



(17) 



In the small time limit it goes ballistically as t 2 as expected, while the mean displacement 
(x(t)) is zero (no atomic current on an un- accelerated lattice). 



2.2. Accelerated lattice in an external AC drive 

When an acceleration is imparted to the lattice, the Block bands are broken up into 
Wannier-Stark Ladder (WSL) of states with level spacing equal to Fd, in other words 
lattice Hamiltonian has a systematic bias. There is a constant energy mismatch between 
the successive site energies (Figure 3) . Now consider that this system of discrete energy 
states is present in an external AC laser field. We describe this physical picture by a 
tight-binding one-band Hamiltonian 

= -E cosu;tJ2[\l)(l + 1| + \l + +^2al\l)(l\, (18) 
1 1 

The wavelength of external standing wave AC field is assumed to be much longer than 

lattice period, so as to subtend a relatively strong dipolar matrix-element between the 
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Figure 3. System with discrete energy levels (WS Ladder) with level spacing Fd 

neighboring WS states of the ladder. So, in the Hamiltonian, the term E cos ut (time 
dependent drive of amplitude E and circular frequency uj) acts as a nearest-neighbour 
transfer matrix element. It may be noted, that in the limit uj = 0, this simulates the 
usual transfer matrix element —V = —E . As before, the time evolution of the reduced 
density matrix (in matrix elements) is given by, 

dp mn _ iE _ _ _ _ .a _ 

^ ^_ [Pm+l,n ~r Pm—l,n Pm,n—1 Pm,n+l\ % ^_ y^Pmn ^Pmn) 

~ Tt 1 - Smn]Pmn (19) 

The quantities ^, and 7 have a dimension of time -1 . We define t = 

<5(dimensionless acceleration) = Jr-, r(dimensionless decohering factor)= to'y, and r 
(dimensionless time)= ^. With this we have 

dpmn • r 1 

d r = 2 c OS^r[p m+ i jn + p m -l, n - Pm,n-1 - Pm,n+l\ 

i5(jnp rnn fip mn ) T[l ^mn]/'rmi (20) 

After applying the rotating wave approximation with p mn = p mn e - lS ( m - n ) T ; the 
evolution of reduced density matrix is, 

dpmn _ i r i(M)rr- _- 1 

q ~ 2 [Pm+l,n l J m,n-H 

+ e" !(8 ^ )T [pm-l,n - Pm,n+l]] ~ Tp mn [l - 5 mn ]. (21) 

Here, we have defined 9 = ujt , r = t/t , ( t = h/E ), and A = 9 — 5 as the detuning 
between drive frequency uj and scaled energy level spacing 5. Note that 

p(h,k 2 ,T) = J2pmn(T)e-^e mk2 

m,n 

p m ,n+i -> e~ tk2 p(k 1: fc 2 , r) , p m+ i, n -> e lkl ~p{k x , k 2 , r) 
Pm,n-i -»■ e lk2 p(h, fc 2 , r) , p m -i, n -> e~ tkl p(h, k 2 , r). 



(22) 
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Figure 4. Plot of mean-squared displacement (x 2 (t, T)) vs scaled time r and 
dimensionless decohering factor T = 7oj^i (proportional to inverse of radiative 
lifetime of excited state). One can clearly identify 'a peculiar transition point' near 
T = 0.5, below the transition point (x 2 (t, T)) increases with increasing T and above the 
transition point mean-squared displacement decreases with increasing T. Oscillations 
smooth away as T increases. Here A = 0.5 



p(ki,k 2 ,r) + 



2tt 



p(k x -q,k 2 - q,r)dq. 



Performing similar co-ordinate transformations p 
defining p(ki, k 2 , r) = g(p, u, r), we have 
dg(p, u, t) 



(h + k 2 )/2 , u 



(23) 

k 2 — ki and 



+ 



dr 
2^ 



= [2i sin(p + At) sin(u/2) — T] x g(p, u, r) 
Q(p ~ Qi u -> i~)dq. 



(24) 



The solution of first order P.D.E. (equation 24) is given in appendix B, we finally obtain 



(At)) 



r 



r 2 + a 2 



A 2 -r 2 
(A 2 + r 2 ) 2 



e rr cosAr} 



2rA 



(A 2 + r 2 ) 2 



e rr sinAr. 



(25) 



Noting that, T (dimensionless decohering factor) = 7-^, r(scaled time)= t^, and 



A (dimensionless detuning)= ^F^- The above equation (equation (25))is an important 
result of the present work. The mean squared displacement from above equation is 
plotted in Figures 4 to 7. For two special cases of interest in long time limit, equation 
(25) gives: 

(A) On-resonance, i.e., (detuning) A = = 0, 



(At)) 



El 
ft 2 7 



t , (diffusive) 



(B) Off-resonance (finite detuning) A 7^ 0, 



(At)) = 



Eh 



h 2/ y 2 + (hu — a) 2 
which indicates diffusion, but with a diffusion constant 

Eh 



t , (diffusive — controllable) 



D 



El 



2[h 2 1 2 + (hw - a) 2 ] 2h2 



7o ( 



U! L —U>0 



[(w—uib)(wl-uq)] 2 



70 



(26) 



(27) 



(28) 
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Figure 5. Enhancement of diffusion with increase in T. The top most curve is for 
r = 0.3, central for T = 0.2, and lowest for Y = 0.1 (for Y < 0.5). With constant 
detuning parameter A = 2. As Y decreases, the oscillations in the mean-squared 
displacement increases, but after a long time oscillations vanish and the mean-squared 
displacement goes linearly with time as it should. Oscillations are a signature of WSL 
due repeated reflections from nearby states. 



<x 2 (r) > 




Figure 6. Suppression of diffusion due to increase in Y (for Y > 0.5). The top most 
curve is for lowest damping constant Y = 0.5, central for Y = 1, and lowest for V = 2. 
With constant detuning parameter A = 0.5. 



tunable with the external derive of frequency u. This is one of the main results of 
this work. The energy-level spacing a = Fd = ush between the WSL states can 
be controlled by the imparted acceleration for optical lattice case and by electrostatic 
field E as a = eE.a for semiconductor superlattice case. Thus a and u acts as 
control parameters in an experiment. Diffusion coefficient becomes maximum in the 
on-resonance case: 

oj c = j: = = u B (Bloch frequency) 



with 



D P 2 El 
2 7 2 7o £ 2 



n 



(29) 



Which is exactly the equation (17), identifying E with V. 




Figure 7. The effect of detuning A ~ Eg a on mean-squared displacement. The top 
most curve is for resonance case, no detuning A = 0, central for A = 0.5, and lowest 
for A = 1. With constant T — 0.08. As the detuning goes up, the oscillations in 
the mean-squared displacement increases, but we have the same expected evolution, 
short time, r 2 rise, then oscillations, and after a long time oscillations vanish and the 
mean-squared displacement goes linearly with time. 



3. Discussion 



We have studied some transport properties of cold atoms in an accelerated and 
harmonically driven optical lattices in the presence of decohering effects due to 
spontaneous emission. The novelty of the work is a tunable diffusion coefficient and 
the possibility to control the diffusve transport by external control parameters. We 
consider a practically important case, in which our system (an atom on an tilted optical 
lattice) is present in an external AC drive with wavelength longer than lattice period, so 
as to subtend a strong dipolar matrix element between the neighboring WSL states. We 
have obtained several interesting results about mean-squared displacement (Figures 4 
to 7). Our main result is a tunable diffusion coefficient (equation (28)), which becomes 
maximum when uac = ^Bioch, (equation (29)), (i.e., the drive frequency is equal to 
Bloch frequency). Similar effects has been observed for ultra cold Bosonic atoms, i.e., a 
system of ultracold bosonic atoms in a tilted optical lattice can become superffuid (Mott- 
insulator to superffuid transition) in response to resonant AC forcing [221 E31 EH 125] . 
Underlying mechanism is the dynamical disappearance of the energy gap due to static 
tilt and external AC drive. There is a possibility to study the effect of decoherence on 
this kind of AC indued flow. 

This diffusion maximization of a single atom should be contrasted with 'localization' 
(dynamical localization of dynamical disorder [26| 127] and Anderson localization of static 
disorder [28J). Dynamical localization is related to the suppression of electron transport 
in a Bloch band driven by an AC electric field, which is explained on the basis of Floquet 
theory of Bloch band collapse [29] . Here we report an opposite effect, a kind of 'maximum 
non-localization'. The physical mechanism under action appears to be the simultaneous 
effects of decoherence (due to spontaneous emission) and dynamical disappearance of 
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energy gap between nearby states of WSL, when drive photon energy is equal to the 
WSL gap energy Because these are the two main physical effects under action. 

For numerical estimation, consider that we have Na atoms (consider very dilute 
atomic sample so that one can neglect atom-atom interactions) on an optical lattice 
created by using two counter propagating laser beams of wave length A = 852nm[25j. 
The resulting lattice spacing is d^ = A/2 = 0.426/im. It is accelerated at a rate 
a = di-^Au = 1000m/sec 2 . Here Au is the difference between the frequencies of two 
lattice beams. The diffusion will be maximum for cuac — ^Bioch = — ma ^ L — 154kHz, 
and the separation (a = Fd) between nearest WSL states will be of the order of lOOpeV. 
The depth Vq of the resulting periodic potential can be QE rec , E rec = ^. . For the 
estimation of diffusion coefficient at resonance between WSL and external AC field, but 
in the far detuned regime between atom and optical-lattice laser field, where only virtual 
transitions can occur, and one can ignore the internal structure of the atom(adiabatic 
elimination procedure), we take Rabi frequency Q = 2pE^jfi = 2tt x 34.5 x 10 7 Hz, 
detuning 5l = 2-k x 5.4 x 10 9 Hz, rate of spontaneous emission (for sodium atoms) 70 = 
6.2 x 10 7 Hz[18\. With this D = ^[^f] 2 = 1.976 x 1CT 6 Sec (in terms of scaled units 
d = l(unit lattice spacing), Eo(amplitude of transfer matrix element) = 1, h = 1). 
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Appendix A: No acceleration and no AC Drive 

We take the time(scaled) laplace transform p(p,u,s) = / °° e~ ST p(p,u,r)dr of equation 
(14). and get 

sp(p, u, s) - p(p, u, t = 0) = (j)(p, u)p(p, u, s) + Tx(u, s). (30) 
We want to calculate the value of p(p, u, t = 0). We know that 
p m ,n(t = 0) = (C* m (t = 0)C n (t = 0)) = 5 mfi 5 n , 



p(h,k 2 ,t = 0) = J^PmAt = 0)e- mk ^e ink2 = J2 6 m,o5 n ,oe- tmkl e mk2 

m,n m,n 

p(p,u,t = 0) = J2<Q<v e ~ im{P+U/2)em{P ~ U/2) = L > 

771,71 

(31) 

with this we get 

p{p, u, s) = — -. 32 

s - <p{p, u) 
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Summing the above equation i.e., equation (32) over p. 



Tp(vu s) = T 1 + mU,s) 



1 

27 



p(p,u,s)dp 



[l + Tx(u,s)} 
2vr 



/7T 
-7T 



dp 



s - (f>(p,u)' 



By re- arrangements we get 

7 



i r 



dp 



1 - IV ' 2tt s - u) ' 
Now, we want to find the mean and mean-squared displacement, as we know 

i r ~ 

X(u,s) = — J p(q,u,s)dq = 

J> m ,„(s)e^ m+n )V(™-- >*dp 

m,ri 



(33) 



(34) 



1 

27 



1 ^ 

" 27 



e i(n ' m)p dp. 



(35) 



The mean displacement in s-domain is given by Yin n Pn,n(s), noting that 
dx(u, s) 



We have 



du 



(x{s)) = i 



-i^2np n<n (s)e~ 

n 

dx(u,s) 



du 



u=0 



Similarly, mean-squared displacement is given by 

~d 2 x{u,s) 



(x 2 (s)} = ^Vp n , n (s) = 



du 2 



(36) 



(37) 



u=0 



Now differentiating equation (34), w. r. t. u and finding the differentials of the integral 
I at u — 0, using equation (36) and equation (37), we obtain 

(x(s)) = , 

(*(*)> = 0. (38) 
1 



After inversion, we get the mean-squared displacement (equation (15)) 



(39) 
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9g( ^ ,T) = [2i sin(p + At) sin(u/2) - T] x g(p, u, r) 

r r 

+ 2^J ^P~^ u ^) d( l- ( 40 ) 
The solution of first order P.D.E (equation (40)) is 

g(p, u, r) = Te-^ U ' T) J e^ (p '"' T) x(«, r)dr 

+ C 1 ( P ,u)e~^ p ' u ' T \ (41) 
In the above, we have defined 

i r i r 

X(u,r) = — J g(p- q,u,r)dq = — J g(q,u,r)dq, (42) 

ip(p, u,t) — — J [2i sin(p + Ar) sin(u/2) — F]dr 
2i 

= — cos(p + Ar) sin(w/2) + Ft. (43) 
Summing over p, 

g(p, u, r) = T e-^ p ' u ' T) I e^ T) x(u,r)dT 
p v 

+ J2 c ^p> u ) e ~^ p,u,T) - ( 44 ) 



^ ^ dpg(p, «, r) = £ jT dpe-^) 



x 



J e lp(p ^ T) X (u,T)dT+ J dpC 1 (p,u)e-^ p ' u ' T \ 



(45) 



X(u,r) = ^-J\^ p ^hdp 

+ 27 \* C ^u)e-^ p ^dp,h = j e^ p ^x(u,r)dr. 

(46) 

To calculate Ci(p, w), we put r = in equation (46), and use the initial condition 
equation (7), i.e.,p mn (t = 0) = p mn (t = 0) = 5 m0 5„o, we have 

0(P,U,T = 0) = Y,Prnn^)e- lm{p - Ul2) e miP+U/2) 

m,n 

= J2 6 ^ 6 noe~ im(p ' u/2) e in(p+u/2) = 1. (47) 



x(u,T)dr 



Controllable diffusion of cold atoms 
So, 

Ci(p,«) = e^> u >°) -TI 4t ,I 4t = 
Equations (46), and (48) gives 

'' u ' T) dp, h = J e^'"' r) x(w, r)dr, 
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(48) 



T=0 



It 



x e 



-<p(p,' 



I At — 



J e^ p ' u ' T) X(u,T)dr 



T=0 



2i 



(p(p, u,t) = — cos(p + At) sin(w/2) + IY. 



(49) 



As we have, (x(r)) 



; 9x{u,t) 
du 



In order to calculate the mean displacement, 



we differentiate the above integral equation for reduced-transformed density matrix 
equation (49) w. r. t. u and set u — 0. Noting that 

i r 

X{u,T) = — J g(q,u,r)dq 

= — [ V Pmn {r)e iS{m - n ^e- im{p - ul2) e m{p+u/2) dp, 
X(0,r) = i- f J2pmn(ry 5{m - n)T e« n - m ^dp 

~ 7r m,n 

E6~ o (r) e i5( . m - n > = I 
u mnl J mn \ ' y ± - 



We finally obtain 

(x(r)> 



du 



0. 



(50) 



(51) 



u=0 



To calculate mean-squared displacement = — 9 



u=0 



we solve equation (49) by 



doubly differentiating it w. r .t. u and then setting u — to get 



'd 2 x{u,T)~ 


= r e - rr , 


/> 


"9 2 x(m,t)" 


du 2 


M=0 ^ 




du 2 



(A 2 + r 2 ) 

-2tt 



(1 - e~ ir cos At) + 



M=0 

e~ 





"9 2 x(m,t)" 


dr 




du 2 


u=0 



A(A 2 + r 2 ) 

T = 



Tt sin At 



(52) 



Equation (52) is solved by Laplace Transform method. In equation (52), we define 



du 2 



u=0 



= /(r). With this, f(r) takes the following form 



f(r) 



f(r)dr 



+ 



cos At 



+ 



sin At 



A 2 + r 2 A 2 + r 2 A(A 2 + r 2 ) 



+ constant. (53) 



Controllable diffusion of cold atoms 



15 



Differentiating the above integral equation, we get 

d HA=Vf(r\ ^ TC ° SAt AslnAr (*a\ 

dr J{) A 2 + r 2 A 2 + r 2 (A 2 + r 2 )' 1 ' 

With the initial condition f{r = 0) = 0, i.e., mean-squared displacement is zero at time 
t = 0, the above equation can be readily solved to get the mean-squared displacement 
as given in equation (25). 
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